Hopping Transport in the Presence of Site Energy Disorder: 
Temperature and Concentration Scaling of Conductivity Spectra 
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Recent measurements on ion conducting glasses have revealed that conductivity spectra for 
various temperatures T and ionic concentrations n can be superimposed onto a common master 
curve by an appropriate rescaling of the conductivity and frequency. In order to understand the 
origin of the observed scaling behavior, we investigate by Monte Carlo simulations the diffusion of 
particles in a lattice with site energy disorder for a wide range of both T and n values. While the 
model can account for the changes in ionic activation energies upon changing n, it in general yields 
conductivity spectra that exhibit no scaling behavior. However, for typical n and sufficiently low T 
values, a fairly good data collapse is obtained analogous to that found in experiment. 
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I. INTRODUCTION 

Electrical conduction processes in disordered solids 
show a remarkable similiarity with xespect to their low- 
frequency dielectric loss behavior. El For frequencies uj 
smaller than some crossover frequency t~^, the real 
part cr'(w) of the complex conductivity a(uj) = cr'{uj) + 
ia"{uj) is constant (cr'(w) = cr(0) = tTdc for ujt <^ 1), 
while for higher frequencies ujt > 1, (7'{lu) increases 
monotonously with w, until at typical phonon frequen- 
cies above 100 GHz vibrational contributions become 
dominant. The dispersive part can be characterized by 
approximate power laws, cr'(w) ^ w*, with exponents 
s ~ 0.6 — 0.8, where s tends to increase weakly with in- 
creasing frequency (at fixed temperature T) and decreas- 
ing temperature. In fact, at low temperatures or high 
frequencies the exponents s becomes close to one. Alter- 
natively, this general behavior may be discussed in terms 
of the complex dielectric function i{ui) = €'{Lo)+ie"{uj) = 
Coo + icr(ijj) / tQtoo, where eq is the vacuum permittivity 
and the high frequency dielectric constant. 

For ionically conducting glasses in particular, it was 
pointed out by Isardo already in 1962 that r^^ and (t(0) 
are proportional with an almost universal constant of 
proportionality. More precisely, it was found that the 
peak frequency LOp in the "dielectric polarization part" 
[cr'(a;) — cr(0)]/w of the loss spectrum is proportional to 
cr(0). However, since any reasonable definition of r^^ 
yields values comparable to Wp, we do not distinguish be- 
tween ujp and here. A closer inspection of the relation 
between and ct(0) was done later independently by 
Barton, Nakajima, and Namikawaa on a variety of ionic 
glasses with the result 



cr(0)T = peoAe' , 



(1) 



where p is a material dependent constant of order one 
and Ae' = e'(0) — too = — £(7^^ l™w^o o'"(a;)/w is the di- 
electric loss strength. Equation (|^) is commonly referred 
to as the BNN relation. 

The question if and to what extent these common fea- 
tures of conductivity spectra in disordered solids can be 
cast into a universal scaling form has been debated for a 
long time. In fact, if r is the only relevant time scale in 
the system and Ae' oc n/T with n being the charge car- 
rier concentration (this is the simplest behavior one can 
expect based on the fluctuation-dissipation theorem) , the 
BNN relation suggests a scaling form 



a(0) 



F{ujt) 



aiO)T 



(2) 



where a is a constant of order unity. This type of scal- 
ing behavior will henceforth be referred to as the BNN 
scaling. 

Various theories of charge carrier transport in dis- 
ordered systems have been developed over the past 
30 years, and already within the classical pair 
approximationP i^nd the continuos time random walk 
(CTRW) modeltl the possibility of scaling was consid- 
ered. The problem has been investigated in more detail 
then by a number of authors. Summerfield suggested,El 
based on the extended pair approximation (EPA) ,13 that 
the real part cr'{uj) in doped crystalline and amorphous 
semiconductors should follow eq. (^) (but without consid- 
ering the dependence of r on the charge carrier concen- 
tration n). Tests of this scaling ansatz for data on evap- 
orated amorphous germanium (for 51 K < T < 102 K) 
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and doped crystalline silicon (for 1 K < T < 10 K) in 
a frequency range 10"'^ — 10^ kHz, showed a good data 
collapse for various temperatures and indeed almost the 
same scaling function i^(ti>T)|-for doued crystalline silicon 
and amorphous germanium.Ei HuntEi suggested a scaling 
behavior for the real part cr'{Lo) using an "augmented pair 
approximation" at intermediate frequencies and a "clus- 
ter theory" at low frequencies, and discussed the degree 
of univpEsality of the BNN relation. 

Dyrellj studied a special electric circuit model resem- 
bling a discretized form of Maxwell's equations for a 
medium with spatially fluctuating dielectric function. By 
treating this model in the effective medium approxima- 
tion (EMA), he showed that at low temperatures the 
complex conductivity o'(w) can be written as 



a{0) 



log 



a{0) 



(3) 



where r oc 1/Ta{0). This low tempetature limit of 
the EMA was first derived by BryksirO and verified 
to hold true for various distributions of the impedances 
in the electric circuit model,t2l thus providing a possi- 
ble explanation for the experimentally observed "quasi- 
universality" . Bryksin's equation (|^) also results f"mi 
an EMA treatment of a microscopic hopping modetj'tj 
with uniformly distributed ^fandom free energy barri- 
ers at low T (see also Rcf.til). This treatment yields 
a crossover time t oc n/Ta{0), in accordance with the 
BNN scaling (H). Note that Bryksin's equation pre- 
dicts a unique scaling function F(z) being determined 
by F{z) log F{z) = iz. 

The Bryksin scaling function was shown to describe 
the normalizjcd real part a'{u!)/cr(0) of various materials 
fairly well,E3 both for electronic conductors (as e.g. n- 
doped crystalline silicon or amorphous germanium) and 
ionic conductors (as e.g. sodium silicate glass or the 
glass former 0.4CaNO3 • O.6KNO3). It must be noted, 
however, that the EMA fails to predict the correct dc- 
conductivity cr(0). The low temperature activation en- 
ergy i?a oc — log(TCT(0)) resulting from the EMA is not 
in accordance with thcr-walue obtained from a critical per- 
colation path analygdSjta which is known to become exact 
at low temperaturesO (see also the results in Sec. Ill be- 



low). Only with a rescaled dc-conductivity may the EMA 
be regarded as a valuable concept to describe conductiv- 
ity spectra in strongly disordered systems. 

For a long time the BNN scaling (||) has never been 
tested with respect to changes in charge carrier con- 
centration n, until Roling et al. recc nthf , studied (7(0;) 
for various ionically conducting glassealJO (see also the 
recent study of Sidebottomcj) . In these materials, in 
contrast to most other systems, the concentration of 
charge carriers (mobile ions) can be varied to a large 
extent. Roling et al. found that the real parts of con- 
ductivity spectra in sodium borate glasses with compo- 
sitions xNa^Q • (1 — a;)B203 (and also other types of 
ionic glasseslla) can be superimposed fairly well onto a 



single curve0 not only for various temperatures, but 
also for various modifier contents x n. However, 
it was found alsal3 that the dielectric strength scales 
as Ae' oc x^^^ n^^^, in contrast to the assumption 
Ae' oc n/T made when deriving the BNN scaling (||). 
This implies that the imaginary part cannot obey the 
BNN scaling because of causality requirements (Kramers- 
Kronig relations). Moreover, in RefH3 alkali germanate 
glasses have been studied over n values varying by a fac- 
tor of ten, and it was found that the BNN scaling has 
to be modified in order to obtain a data collapse for this 
large concentration range. We therefore have to consider 
the observed data collapse in cr(w) as an effective one. 

The effective scaling behavior is nevertheless quite 
remarkable, since the activation energy i?a of the dc- 
conductivity (t(0) (and hence the activation enHcgy for 

with 



oc cr(0)) usually decreases logarithmicall; 



ionic concentration n, = A — B\ogn, for typical 
n « lO'^^ — 10^^ cni^-^. Accordingly, one finds a quite 
complicated explicit dependence of the crossover time t 
on n, T n^~^/^'^^ . 

It is important to note that these experimental findings 
for ionically conducting glasses cannot be explained by 
the random free energy barrier model, since in this model 
the conductivity is proportional to n but does not steeply 
increase with n as (t(0) ~ n^/^^'^ . The perhaps simplest 
approach that allows one to account for a concentration 
dependence of i?a is a model of particles iippping in a 
lattice with site energy disorder (see Refs.EJ and be- 
low), where only one particle can occupy a given lattice 
site (Fermi statistics) . In order to see if such a model also 
yields the BNN scaling (|^) , we perform Monte Carlo sim- 
ulations for Gaussian distributed site energies in a range 
of concentrations n and temperatures T both extending 
over almost two orders of magnitude. We find that the 
model can account for the logarithmic dependence of Eg_ 
on n, but in general exhibits no effective BNN scaling 
behavior corresponding to (|^). However, for typical ionic 
concentrations and sufficiently low temperatures an ef- 
fective BNN scaling behavior is found in agreement with 
experiment. 



II. HOPPING MOTION AND 
TRANSPORT QUANTITIES 

The jump motion of particles in a lattice with site en- 
ergy disorder is a standard modeLfor describing hopping 
transport in disordered systems .EZl We choose here a sim- 
ple cubic lattice with spacing a and assign to each lattice 
site i independently a random energy Si drawn from a 
Gaussian distribution tp{e) with zero mean and variance 
(Tg . The particles jump among nearest neighbor sites and 
the jump rate from site z to a vacant neighboring site j is 
W{£i,£j) = i/min(l, exp[— (£j — £i)/k-QT]) with v being 
an attempt frequency (Metropolis transition rates). If 
the neighboring site j is occupied, the jump rate is zero. 
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The size of the lattice is Li^ and periodic boundary condi- 
tions are used. At equihbrium, the probabihty for a site 
with energy e to be occupied is given by the Fermi distri- 
bution /(e) = [l-|-exp[(e — ^)/fcBr]^^, where the concen- 
tration c = na^ per lattice site determines the chemical 
potential through c = {a/ L^Yli fi^i) — J de^p{e)f{e). 
For T 0, /(e) approaches the step function /(c) = 
9{ef{c) — e) {0{x) = 1 for x > and zero else) with the 
Fermi energy ef (c) given by 

erf(ef(c)/2i/V) = 2c-l, (4) 

where erf(-) denotes the error function. 

In order to analyse the mobility of the particles for var- 
ious concentrations c and temperatures T we determine 
the mean square displacement (r^(i))|-Cjf a tracer particle 
by means of Monte Carlo simulations Ba and calculate the 
time-dependent diffusion coefficient D{t) = {r^{t))/6t. 
D{t) approaches the short-time diffusion constant Dgt 
for t — > 0, and the long-time diffusion constant D^o for 
t — > oo. The frequency-dependent conductivity (t(w) is, 
to a good approximation,E3 given by 

2 2 roo 

'7(t^) = -^^ lim / (r2(t)) exp(it^t-ei) dt, 
fcsT 6 e^O+/o 

(5) 

and approaches cr(0) — f^Doo and a{oo) = f^Dgt in 
the low and high frequency limits, respectively. 

The short-time or high-frequency dynamics is governed 
by the mean jump rate of a particle, 

(^W) = j de MM J de' ^{e') [1 - f{e')] W{e, e') , 

(6) 

which determines Dgt = (W)a^/6. For fcsT ^ |ef(c)| we 
can use the high temperature expansion of the chemical 
potential, ^ = A:BT[log[c/(l - c)] (c - l/2){a,/kBTy + 
0{a^/kBT)'^) to evaluate eq. (^ and obtain 

At low temperatures /sbT ^ |£f(c)L we can replace 
/(c) by the step function in eq. (^) and it follows 
that {W) = {Av/c) exp[(cre/fcBT)2] evicE- evicE+ with 
E± = 2-i/2(a,/fcBr±Cf(c)/ae) (erfc(-) - I - erf (•) de- 
notes the complementary error function). For k-oT ^ 
min(crg/|ef (c)|, |cf(c)|) in particular, we then find, by us- 
ing erfc(2:) ^ exp(— a;^)/7r^/^a; for x +oo, 

W-^(^)'exp[-ef(c)V-3. (8) 

Note that (W) is not thermally activated in this low 
temperature limit. This is a desirable feature of the 



model, since the high frequency conductivity is indeed 
only weakly temperature dependent in ionically conduct- 
ing glasses (and other disordered hopping systems). 

By contrast, the long-time diffusion constant or dc- 
conductivity is thermally activated at low temperatues, 
Doo oc cr(0) oc exp(— iJa/fcer). The activation energy 
i?a can, at sufficiently low c (['0(ef (c))fc^]^/'^ ^ 1), be 
calculated from the critical path analysial3 as the differ- 
ence between the critical energy Cc and the Fermi energy 
Cf (c), i.e. -Ea = £c — ef(c). Here, Cc is the highest energy, 
which at least has to be encountered by a particle when 
it wants to move between two (far distant) sites close , 
to the Fermi level. According to percolation theoryEy'Ell 
fj^ V(e) de = [l+erf(ec)]/2 = pc, where Pc = 0.3117 de- 
notes the percolation threshold in the sc-lattice. Solving 
numerically for Cc we obtain Cc = — 0.491cr£ and hence 

£;a =^ -0.491ae-ef(c). (9) 

For c 0, £f{c) — oo, such that the 

implicit equation erf(ef (c)/2^/^(T£) = 2c — 1 for 
ef(c) can be simphfied to exp(— [cf (c)/(2^/^(Te)]^) ~ 
-27ri/2^[ef(c)/(2i/V£)], which yields ef{c)/cr, ~ 
-2i/2[- log(27ri/2c[- log(2^i/2c)]i/2)]i/2. 



III. RESULTS AND DISCUSSION 

Figure 1^ shows the Monte Carlo results for D{t) (in 
units of Dq = va?) and the real and imaginary parts 
<j'(u)) and cf"{uj) of the conductivity (in units of cro = 
ne^D[j/k-QT) being calculated from D{t) by numerical a 
Laplace transformation according to eq. (|^). In the up- 
per figures (Figs. ||a-c), the curves refer to a fixed con- 
centration c = 0.02 and various (inverse) temperatures 
ranging from a^/k-QT = 2 to 7, while in the lower figures 
(Figs. Qd-f) the temperature is held fixed, ae/k^T = 6, 
and the concentration is varied between c = 0.005 and 
c = 0.16. The overall behavior of the complex con- 
ductivity (t{uj) compares fairly well with that found in 
the experiments. There exists a crossover time r, such 
that the real part cr'i'jj) is almost constant for lot <^ 1, 
<j'{io) ~ c(0). For LOT > 1, (7'{uj) increases monotonously 
with Lo until for uo ~ Whf it reaches the high frequency 
plateau being determined by (W) (see Sec. |l|). In real 
materials this high frequency plateau is usually masked 
by phonon excitations, but may bejiesolved by substract- 
ing the vibrational contributions .eHI The dispersive part 
of the spectrum for r^^ ^ ^ w^f can be described by 
an approximate power law at sufficiently low tempera- 
tures {ae/ksT ^ 4, see Fig. |l|b) with typical exponents s 
in the range 0.5-0.8. The imaginary part (j"{uj) behaves 
similar to <j'{lo) in the dispersive regime t^^ <C t^J <C com, 
i.e. (t"{lo) ^ Lo^ , while in the dc- and high frequency 
limits it rapidly approaches zero, since the phase shift 
between current and electric field approaches rapidly 
and TT, respectively. 
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FIG. 1. Plot of (a),(d) D{t)/Do vs ut, (b),(e) (T'(aj)/cro vs uj/v, and (c),(f) cr"(aj)/(To vs uj/u (ao = ne^Do/kBT). In (a-c) the 
concentration is fixed, c = 0.02, and the temperature is varied, a^/k^T = 2 (full squares), 3 (open diamonds), 4 (full triangle), 
5 (open stars), 6 (full circles), and 7 (open squares), while in (d-f) the temperature is fixed, a^/k^T = 6, and the concentration 
is varied, c = 0.005 (full stars), 0.01 (open circles), 0.02 (full squares), 0.04 (open diamonds), 0.08 (full triangles), and 0.16 
(open stars). The solid lines in (a) and (d) represent the fits used for the Laplace transform. 



Figure ga displays Arrhenius plots of a^cT oc cDoc/Dq 
for various c. The linear behavior of the curves at low 
k-eT <Si airi^ee the dashed lines in the figure) allows us to 
determined the activation energy E^^ whose dependence 
on c is shown in Fig. |^ (in a semi- logarithmic plot). This 
dependence is in good agreement with the behavior pre- 
dicted by the critical path analysis (see the solid line in 
Fig. ||b). As indicated by the dashed line in Fig. ||b, the 
data for c ^ 0.01 can be approximated by a logarithmic 
dependence E^ja^ = A — Blogc, in accordance with the 
experimental findings. We also see that for low concen- 
trations c < 0.01, the data start to fall off the straight 
dashed line. This is expected to occur in view of the low-c 
limit of £f(c) discussed after eq. (^. The behavior is not 
in contradiction with the experimental situation, how- 
ever, since in ionically conducting glasses the logarithmic 
c-dependence of i?a is typically found at relatively high 
ionic concentrations (corresponding to network modifier 
contents of 5-40 mol %, which arc most interesting for 
applications). 

Let us now discuss as to what extent the data can be 
superimposed by a BNN scaling ansatz as in eq. (^. To 
this end the data in Fig. |l| are plotted in a scaled form 
in Fig. ^. As can be seen from Figs. |^a-c, the data for 
fixed c and various T do not collapse onto one master 
curve. The BNN scaling with respect to various concen- 
trations at fixed (low) temperature T is significantly bet- 
ter (Figs. ||i-f), but it is important to note that this ap- 
proximate scaling behavior ceases to be present at higher 




FIG. 2. (a) Plot of cDoo/Dq vs a^/ksT for concentrations 
c — 0.001 (open triangles) 0.005 (full stars), 0.01 (open cir- 
cles), 0.02 (full squares), 0.04 (open diamonds), 0.08 (full 
triangles), and 0.16 (open stars). The, solid line indicates 
the high-temperature activation energyu Eo = a^/n^^^. The 
dashed lines are the fits used to determine the activation en- 
ergy Ea,. (b) Plot of the activation energy E^ja^ vs c. The 
open circles represent the values obtained from the Monte 
Carlo simulation shown in (a), while the solid line marks the 
result from the critical path analysis. The dashed lined is 
a fit with respect to a logarithmic dependence of E^ on c, 
E^ = A-B log(c) with A 0.23 and B 0.37. 
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FIG. 3. Scaling plot of (a),(d) D{t)/Doo vs t/r, (b),(e) a'{Lj)/adc vs ujt, and (c),(f) a" [lu) / a^^ vs lot, for the data presented 



Fig. 0. 

T, a^/ksT < 5. We conclude that the hard core lattice 
gas with Gaussian site energy disorder in general exhibits 
no BNN scaling. 

However, the data collapse for the lowest temperatures 
{a^/k^T = 6 and 7) in Fig. ^a-c is fairly good. In or- 
der to see if this holds true also for other c values, let 
us first estimate what is a typical range of ionic concen- 
trations used in the experiments. For this purpose one 
might consider sodium silicate glasses a;Na20 • (1 — a;)Si02 
with modifier contents x between 10% and 40%. By as- 
suming about 10 sites per silicon to exist, we would have 
c = 2a;/[10(l - x)] ^ 0.01 and c ^ 0.13 for x = 0.1 and 



X = 0.4, respectively. For these c values we always found 
an approximate BNN scaling for ae/k-gT > 6. A scal- 
ing plot combining various low T and typical c values is 
shown in Fig. ^. The quality of the approximate data 
collapse is not as good as the one found in experiment, 
but might become better at lower T. In fact, when con- 
sidering typical activation energies for ionic conduction 
in glasses of 0.3 to 1 eV, one should assume similar val- 
ues for (Te (see Fig. ^d), and accordingly a^/k^T < 10 
at room temperature (this, however, is true only when 
neglecting the Coulomb interaction between the mobile 
ionf^l). 
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